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At ultracold temperatures, the Paul! exclusion principle suppresses collisions 
between identical fermions. This has motivated the development of atomic 
clocks using fermionic isotopes. However, by probing an optical clock transi- 
tion with thousands of lattice-confined, ultracold fermionic Sr atoms, we have 
observed density-dependent coUisional frequency shifts. These collision effects 
have been measured systematically and are supported by a theoretical descrip- 
tion attributing them to inhomogeneities in the probe excitation process that 
render the atoms distinguishable. This work has also yielded insights for zero- 
ing the clock density shift. 

Quantum statistics play a critical role in shaping interactions between matter. This is ap- 
parent in the markedly different behavior of Bose-Einstein condensates (i, 2) and degenerate 
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Fernii gases of ultracold atoms (5). The quantum statistics of atoms can thus be a key factor in 
the choice of an atomic system for a given experiment. Such is the case for atoms at the heart 
of an atomic clock. Simultaneous interrogation of many atoms is favorable for achieving high 
measurement precision. However, when atoms interact with each other, their internal energy 
states can be perturbed, leading to frequency shifts of the clock transition (4, 5). The use of 
identical fermions was prescribed to allow many atoms to strengthen the signal without such 
density-dependent collision shifts (6). Previous experiments seemed to confirm this fact for 
both single-component (7) and two-component fermion mixtures (S). 

However, by probing an optical clock transition with thousands of fermionic Sr atoms con- 
fined in a one-dimensional optical lattice, we clearly observe density-related frequency shifts at 
a fractional precision of 1 x 10~^^. When the light-atom interaction introduces a small degree of 
inhomogeneous excitation, previously indistinguishable fermions become sUghtly distinguish- 
able. This effect causes a time-dependent variation of the two-particle correlation function, giv- 
ing rise to an apparent mean-field energy. The resulting collision effects have been measured 
systematically as a function of temperature, excitation probability, and interaction inhomogene- 
ity. These observations are supported by a theoretical description of fermionic interactions 
including the effect of the measurement process. 

The latest generation of optical atomic clocks such as those based on the ^Sq— ^Po transi- 
tion in fermionic ^^Sr currently offers the highest measurement precision, useful for measuring 
possible atomic interactions (9, 10). In an ultracold dilute gas with a mean field energy, a 
narrow clock transition will experience a density-dependent frequency shift (11, 12) given by 
/lAi/ = {A'Kh^G^'^'^pa)/m. Here m is the atomic mass, p is the density of the atomic sam- 
ple, a is the s-wave scattering length characterizing the atomic interaction, and h = 27rh is 
Planck's constant. G*^^) is the two-atom correlation function at zero distance, which summa- 
rizes the quantum statistics of colliding bodies. For example, G^^^ = for identical fermions 
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and G^^) = 2 for identical bosons in a thermal gas. The Fermi suppression arises from the 
Pauli exclusion principle, which prohibits even-partial-wave collisions between indistinguish- 
able fermions. At ultracold temperatures partial waves higher than s-wave are frozen out (13). 
For atoms excited in our two-level clock system, three possible s-wave interactions exist: those 
between two ^Sq ground-state (\g)) atoms, those between two ^Pq excited-state (|e)) atoms, and 
those between a \g) atom and a |e) atom. Including all possible interactions, the colUsional 
frequency shift at ultracold temperatures is given by Eq. 1 (8, 12, 11): 

Ai/ge = —{GfJageiPg - Pe) + GfJueePe - Cf^aggPg), (1) 

lib 

where Oj^ is the s-wave scattering length for colUsions between atoms in state i and j, and pi 
is the density of atoms in state i. Since indistinguishable fermions do not collide, Ggg = Gee 
= 0. Fermions in different internal states are distinguishable, and for a completely incoherent 
mixture of the two states, G^J — 1. However, if the two-state mixture is prepared by a uniform, 

(2) 

coherent excitation of ground-state atoms, then the fermions evolve indistinguishably and Gge 
= (8). In this case, Auge = 0. 

Two possibilities exist for Ai/^g to deviate from zero. First, the p-wave contribution may not 
be negligible. However, for ultracold atoms confined in a well-characterized optical trap, we 
show experimental evidences and theoretical calculations that conclude that p-wave collisions 
make no noticeable contribution to the observed clock frequency shift. Second, it is impera- 
tive to consider the entire interaction, including the measurement process, when exploring the 
question of whether fermions collide. In fact, the measurement process, such as probing a clock 
transition, may strongly influence the time-dependent variations of G^^\ We show here that an 
inhomogeneous interaction between light and atoms leads to the loss of indistinguishabiUty of 
the fermions, thus making < G*^^^ < 1. 

Although a uniform, coherent excitation of identical fermions maintains G^^^ — 0, and no 
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s-wave collisions occur, if a small non uniformity in the excitation process arises, the atoms 
are no longer completely identical, and G*^^^ > 0. The value of G^"^^ will depend on the degree 
of excitation inhomogeneity. This measurement-induced dynamic variation of quantum statis- 
tics leads directly to a change of the mean-field energy within the ultracold gas, resulting in a 
nonzero Az/^g- It is interesting to contrast the present work with previous results observed using 
an ultracold gas of fermionic ^Li, where the insensitivity of a radio-frequency transition to coUi- 
sional shifts was demonstrated (7,8). It was shown that the fermionic insensitivity to colUsional 
shifts was maintained even when a pure superposition state of the two level system had deco- 
hered. This decoherence would have allowed interactions, but when a uniform rf probing field 
reintroduced coherence to the atoms in a homogeneous manner, (7(2) again became zero, giving 
no coUisional shifts within the measurement precision (14). ^From the current experiment it is 
clear that any non-identical evolutions during the interrogation process lead to the breakdown 
of Fermi suppression; this experiment is sensitive to very small inhomogeneities because of the 
high measurement precision. 

An intuitive understanding emerges from considering two sample atoms in a pseudo spin- 
1/2 system with ground \g) and excited |e) states. Before applying the spectroscopy pulse, the 
atomic system is in a pure, polarized spin state with {i/ji) = 1^2) = \g)- The effect of the pulse 
is to perform a rotation on the Bloch sphere, as shown in the inset of Fig. IB. For a coherent, 
homogeneous excitation, the wavefunction of the system becomes a coherent superposition 
iV'i) = IV'2) = 0!\g) + The wavefunctions of both atoms are identical, = 0, and 
collisions cannot occur. An inhomogenous spectroscopic excitation, such as that caused by 
varying Rabi frequencies for different atoms, results in slightly different rotations on the Bloch 
sphere for the two atoms (Fig. IB inset). Hence we have = a\g) + (3\e) and 1^2) = 

(2) (2) 

j\g) + S\e). The fermions are distinguishable and < <!• The value of depends 
on the amount of inhomogeneity, and its time variation can be explicitly calculated from the 
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anti-symmetrized overlap of the two wavefunctions (details in supporting text (15)): 



Gg\a{t),P{t),j{t),S{t)) = 1 - \a{t)j*{t)+P{t)S*{t)\'. 



(2) 



The resulting collision shift from Eq. 1 is then 



Aiy{t) 



2ha. 
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Before proceeding with experimental results, we first summarize the system under study 
(75). In the ^^Sr optical clock, atoms are trapped in a one-dimensional optical standing-wave 
potential (ID optical lattice). Longitudinally the atoms are confined tightly, with an oscillation 
frequency ~80 kHz. At temperature T = 1 jjK, ~98 % of the atoms occupy the ground- 
state of the trap (n^ = 0.02). The laser probing the clock transition propagates along the lattice 
axis, and spectroscopy is performed in the Lamb-Dicke regime. In the transverse plane the 
confinement is much weaker, with an oscillation frequency — Vy ~450 Hz, and atoms occupy 
a large number of motional states (n^ =^y = 46). Typically, ~2x 10^ atoms are trapped in the 
optical lattice, resulting in 30 atoms per lattice site with a density of 2x 10^^ cm~^ (75). The 
optical lattice is nearly vertically oriented and is operated at the so called "magic wavelength" 
of \l ~813.429 nm (76), where the ac Stark shifts of the ^Sq and ^Pq states are identical. 

With a perfect alignment of the probe laser along the strong confinement axis, assuming 

— * 

cylindrical symmetry, a residual angular spread between the probe and lattice k remains due to 
the finite size of the lattice beam (77). However, an even larger effect occurs if the symmetry is 
broken due to either aberrations in the beam profile or angular misalignment (A6') between the 
lattice and the probe beam. For our trap parameters, we estimate an effective 10 mrad 

(Fig.lA inset). The residual wave-vector projected on the transverse plane leads to slightly 
different excitation Rabi frequencies Vtfi for atoms in different {rix, riy) states (75, 18, 19). For a 
given T, the occupation of a transverse motional state nx,y is given by the normalized Maxwell- 
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Boltzmann distribution. The inhomogeneity in the Rabi frequencies is thus affected by both T 
and A^. 

To calculate the density shift, we return to our two-atom model. Each atom has a slightly 
different Jlfj. For the entire atomic ensemble, we can define an average Rabi frequency Q and 
its RMS spread Af2. To approximate the average density shift, we set l^i = ^2 + Af2 and 
0,2 — ^ — AO for our two-atom model. Thus the time-dependent quantities, a, (3, 7, 5 as 
defined in Eq. 2 are parameterized by Q and AQ {15). At a time t during the spectroscopy 
pulse, the atoms experience an ensemble- averaged shift: 

Az/(t) = ^^gSJ(^+ Al^,^-Al^)(Pfl-Pe). (4) 

This shift evolves during the spectroscopy pulse, and for the final density shift we time average 
Az/(i) over the total pulse length tp- This approximation is valid in the limit that the change 
in Q, due to atomic interactions is much less than Afi. A more rigorous calculation using the 
optical Bloch equations including atomic interactions has also been made. Using our typical 
trap parameters we find the two-atom approximation is valid to within 5 %. The time-dependent 
Rabi-oscillation is only sUghtly affected by atomic interactions, however the effect on the final 
clock shift is obvious. 

For inhomogeneity-induced collision shifts, is important. Atoms in close proximity to 
each other tend to have similar Rabi frequencies, whereas atoms located far apart are more likely 
to experience different excitations (and hence be distinguishable). If tpv^^y 1, the atoms are 
effectively frozen in place, and will experience no density shift. However, if tFi^x,y > 1> atoms 
initially located far apart have time to interact. For the clock experiment requiring high spectral 
resolution, we have = 80 ms and l/z^^.y = 2.2 ms, so collisions will occur. 

To systematically study these effects, we implemented controlled variations of both T and 
A^. To vary T, we perform cooling (heating) of the lattice-confined atoms in three dimensions: 
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Doppler cooling (heating) along the transverse direction and sideband cooling (heating) along 
the longitudinal axis. Simultaneous with the sideband cooling (heating), the atoms are spin- 
polarized by optical pumping in a weak magnetic (B) bias field. Atoms are polarized into either 
the mi;' = +9/2 or mp = -9/2 Zeeman states. The ^Sq — ^Pq clock transition, which is predicted to 
have a natural linewidth of ~ 1 mHz (20, 21, 22), is interrogated using a cavity-stabilized diode 
laser at 698 nm with a linewidth below 1 Hz (25). Spectroscopy is performed in the Lamb- 
Dicke regime and in the resolved sideband limit (24). To ensure that the polarized spin state is 
well resolved from other mp levels, spectroscopy is performed under B ~250 mG, leading to 
a separation of 250 Hz between the = ±9/2 states. A spectroscopy pulse length oftp-^O 
ms results in a Fourier-limited linewidth of ~ 10 Hz. 

After the spectroscopy pulse is applied, atoms remaining in l^f) are counted by measuring 
fluorescence on the strong ^Sq— ^Pi transition. Atoms transferred to |e) are then pumped back 
to \g) via the intermediate (5s6s)^Si states and are also counted. Combining these two mea- 
surements gives us a normalized excitation fraction Pe/iPe + Pg)- The atomic temperature is 
determined using both sideband spectroscopy (25, 75) and time-of-flight analysis. In Fig. lA, 
sample spectra are shown for two different values of T. Once T is measured, the degree of 
inhomogeneity is determined by fitting the decaying Rabi oscillations for the ensemble. In Fig. 
IB, the Rabi oscillation at T = 3 /xK (squares) clearly shows faster dephasing than that of T = 1 
fjK (circles), indicating a larger degree of inhomogeneity. 

Density-dependent frequency shifts of the ^^Sr clock transition are measured using a re- 
motely located calcium optical standard at NIST (9) as a stable frequency reference, which is 
linked to JILA via a phase-coherent fiber network (26). This direct optical frequency measure- 
ment between two optical standards allows fractional measurement precision of a few times 
10~^^ after hundreds of seconds of averaging. To measure the clock center frequency, the spec- 
troscopy pulse is first appUed to atoms optically pumped to the = +9/2 state. In the next 
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cycle, atoms polarized to the m^? = -9/2 state are used. The center frequency is then determined 
by the average of both resonances. The density-dependent frequency shift is determined using 
an interleaved scheme, where the density of the atomic ensemble is varied every 100 s. The 
density is varied by a factor of two. Pairs of such data are then used to measure a frequency 
shift, and many pairs are averaged to decrease the statistical uncertainty. Typically, we lock the 
clock laser near the full- width at half -maximum of each resonance, however the location of the 
lock points is varied to select the desired excitation fraction. 

Spectroscopy is performed using two different experimental procedures. In the first, we 
probe the clock transition from l^f) to |e) (Fig. 2 inset). The intensity of the probe is set to 
produce a 7r-pulse on resonance. This direct scheme could suffer from imperfect polarization of 
the atomic sample, and spectator atoms could be left in other m^;' levels. This scenario could po- 
tentially lead to density-dependent shifts due to collisions between different mp states that are 
not suppressed by the Fermi statistics. The second scheme minimizes this effect by probing |e) 
to \g) (Fig. 2). Here we apply a strong pulse to first transfer the population from l^f) to |e). The 
pulse power broadens the transition in order to decrease the sensitivity of population transfer to 
probe laser frequency, and transfers ~50 % of the population to |e). This first pulse is resonant 
with atoms in one of the mp = ± 9/2 states, hence atoms left in other mp states due to imper- 
fect polarization are not transferred. Subsequently, all atoms remaining in \g) are removed from 
the lattice with a pulse of light resonant with the strong ^Sq— ^Pi transition, without affecting 
the temperature of the atoms in |e). This is confirmed with sideband spectroscopy (15). Finally, 
the clock transition of |e) to \g) is probed with the usual 80 ms 7r-pulse. In both experimen- 
tal procedures, we measure populations in |e) and \g) to determine the normaUzed excitation 
fraction. 

Figure 2 summarizes the measured density-dependent frequency shift as a function of the 
normalized ground-state fraction for two different values of T, 1 /xK (squares) and 3 /xK (cir- 
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cles). The data indicate a clear trend that the density shift decreases under a more homogeneous 
excitation. The solid lines are the expected shifts calculated from the two-atom model. For 
clock operation, it is important to note that near 50 % excitation fraction, for both T, the shift 
goes through zero. 

Of course, as we change T, we vary both the excitation inhomogeneity and the p-wave 
contribution. To estimate the magnitude of p-wave colUsions, we note that the van der Waals 
potential for all three interaction types {gg, ee, or eg) has been theoretically calculated (27, 
21, 28), and the p-wave centrifugal barrier is expected to be greater than 25 fiK. At T ~1 ^K, 
ka <C 1, where k = 2ti/\t. Xt = /i/V27rrnAJBT is the thermal deBrogUe wavelength, and ks 
is the Boltzmann constant. Under these conditions, the ratio of p-wave to s-wave phase shift 
is {bk)%/a, where b is the p-wave scattering length. For gg interactions, the s-wave scattering 
length has been measured (29) for *®Sr, and mass scaling gives agg= 96.2(l)ao for ^^Sr, where 
oo is the Bohr radius. Combined with the van der Waals potential, the p-wave phase shift 
can be determined from the Schrodinger equation. For ^Sq, bgg = -76 Oq, and for T = 1 fxK, 
\{bggk)%gg/agg\ 0.01. Thus, P'Wave collisions for gg are suppressed by over two orders of 
magnitude and are negligibly small. Although the s-wave scattering lengths a^e and a^e have 
not yet been measured and thus cannot directly constrain the values of b^e and b^g, calculations 
based on a theoretical potential predict that these p-wave collisions are similarly suppressed 
relative to s-wave collisions. An exception would be a p-wave shape resonance {13); however, 
this would occur only for a very small range of possible a,,„ and a^e, and the effect would be 
reduced by thermal averaging. We also note that in a trapping potential, k is modified due to 
the zero-point energy of the trap {kzp) and the effective thermal wavevector for collisions is 
given by kr = -|- /c|p)/2. For our trap, kzp ~3.5 fiK, and p-wave collisions are still 

suppressed. The observed density shift scales as age> and for our typical temperatures we 
find values of G^i2 between 0.03 and 0.15, whereas the p-wave scattering length is expected to 
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be ~1 % of age. Hence, inhomogeneity -induced s-wave collisions dominate. In the the unitarity 
limit where k'r\age\ > 1 (a^e is the zero-temperature scattering length), the effective scattering 
length is l/Zcr- For our lattice trap parameters and temperature range of 1-3 //K, this length is 
on the order of -300 oq, which is consistent in sign and magnitude with our observed frequency 

(2) 

shifts, along with the values and uncertainties of G\2 and p. 

To provide further evidence to exclude p-wave contributions, we vary the inhomogeneity 
by misalignment of the spectroscopy probe beam under a fixed T. This also helps rule out 
nx,y,z - dependent residual ac Stark shift of the trap. Typically the probe beam is coaligned 
with the lattice to minimize motional effects. However, by increasing the misalignment (A^), 
we can also increase AQ. Fig. 3 A and B show Rabi oscillations for two different probe beam 
misalignments at T = 1 fjK (triangles and open squares) and 3 fxK (circles and open diamonds), 
respectively. Fig. 3C displays the measured density shift as a function of (Ail/f2) due to probe 
misalignment. For T = 1 fxK, the shift becomes larger with increased AQ/fi. When AQ/Q 
increases further, the 3 /xK data indicate that the density shift becomes smaller. This behavior is 
reproduced by the theoretical curves shown in Fig. 3C, and is illustrated in Fig. 3D. Consider 
two different AQ/Q, both with an average excitation fraction of 0.3. In the first case, for small 
misalignment, we find a spread in the excitation fraction of ±0.2; there is an inhomogeneity 
allowing collisions to occur and we measure a small density shift. In the second case, with 
further misalignment the spread in the excitation fraction increases to ±0.4; there is now a 
larger spread in the Rabi frequencies, and collisions still occur. However, we now have atoms 
with an excitation fraction both above and below 50 % where the shift crosses zero. Hence, 
the collisions of atoms with excitations between 0.3—0.7 will average to zero (this is consistent 
with the density shift going to zero at 50 % excitation, regardless of the inhomogeneity), and 
the final collision shift is due only to atoms with excitation fractions between and 0.3. The 
measured shift for the larger misalignment is therefore smaller. 
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Combining the measurements shown in Figs. 2 and 3 makes it clear that the observed 
density-dependent shifts arise from the change of the quantum statistics G^'^^ caused by the 
inhomogeneous measurement process. Of course, the inhomogeneous effect can be suppressed 
by decreasing the sample temperature and increasing the transverse confinement, or going to 
higher dimension traps. Importantly, for clock operations we have identified that near a 50 % 
excitation fraction the density shift goes to zero. Using these measurements we can now reduce 
the uncertainty of the collision shifts for clock operation (9) to 5xl0~^^. This time dependent 
variation in quantum statistics will also apply to boson-based clocks, where the original G*^^^ = 
2 will decrease to a value between 1 and 2. 
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Figure 1: (A) Sideband excitation spectra for T = 1 /iK (blue circles) and 3 fiK (red squares). 
The spectra are obtained in the resolved sideband limit and have three dominant features, the 
narrow carrier transition and broad red (blue) motional sidebands that are excited when an 
atom is transferred to a lower (higher) motional state during the transition. As the temperature 
of the sample is lowered, the atoms primarily occupy the ground-state and the red sideband is 
suppressed. The temperature of the atomic ensemble can be extracted from a fit of the sidebands 
(25). The inset shows the lattice geometry and excitation scheme. The probe beam and lattice 
are coaligned and copolarized, minimizing the relative spread between ki and kp. However, even 
under the best effort, a small angle A9 between the probe and lattice beams may persist due to 
aberrations and misalignment. (B) Rabi oscillations for temperatures of 1 fiK (blue circles) and 
3 fiK (red squares). For higher temperatures, more motional states are occupied. This leads to a 
larger spread in the Rabi frequencies and faster dephasing of the excitation between atoms. By 
fitting the decay of Rabi oscillations we can determine the degree of excitation inhomogeneity. 
The inset illustrates the dephasing process using rotations on the Bloch sphere. At time a, 
before the excitation, the atoms are in a pure state. At time b, the atoms have undergone two 
oscillations. For the red curve the temperature is hotter and there is a larger spread in Rabi 
frequencies. This is indicated by the increased width of the Bloch vector and dephasing of the 
observed oscillations. At time c, the effect is even more pronounced. 
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Figure 2: Measured density-dependent frequency shift as a function of the final excitation frac- 
tion and temperature. Atoms are initially spin polarized and transferred to ■'Po (|e)) before the 
spectroscopy pulse is applied. The squares (circles) show the measured shift for T = 1 (3) 
fiK. The lines show the calculated shifts using the two-atom model with only a single scaling 
factor. Near ~50 % the shift goes through zero. In the inset, the measured shift is shown for 
atoms excited from the ^Sq {\g)) state for T = 1 fiK (squares), 3 /iK (circles), and 5 fiK (trian- 
gles). However, the magnitude in this case could be influenced by imperfect spin polarizations. 
For both plots, as the temperature is decreased, the inhomogeneity also decreases, leading to 
a smaller collision shift, p is the atomic density of lO^Vcm^. The density dependent shift for 
each excitation fraction is determined using an interleaved scheme where the density is varied 
every 100 s. Pairs of such data are then used to determine the frequency shift. Typical data sets 
include 20 - 30 pairs of density comparison, with the error bars indicating the standard error. 
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Figure 3: Effect of probe misalignment on the density-dependent shift. (A) Rabi oscillations 
are shown for two different values of A6 at T = 1 /xK. The open squares show oscillations when 
the probe is aligned similar to that of Fig. 1 and 2. The solid triangles show a faster dephasing 
when the probe beam misalignment is increased further by 5 mrad. (B) Rabi oscillations for T 
= 3 fiK. The circles show oscillations when the probe beam is aligned similar to that of Fig. 1 
and 2, and the diamonds when the misaligment is increased further by 35 mrad. (C) The density 
shift measured for each misalignment shown in (A) and (B). ^From A9 and T, the spread in 
Rabi frequency AVl is calculated. The lines show the expected shift as a function of A^l for 
T = 1 fxK (solid line) and 3 fiK (dashed line). The inset shows a zoomed-out plot. (D) For 
large misalignments, we observe a smaller density shift. This is described using the rotation 
on the Bloch sphere. As an example, two different values of AD, are shown. On each sphere, 
the average excitation fraction is shown with a solid line, and the spread is indicated by the 
dotted lines. For small misalignments, we have a small spread in Rabi frequencies. As the 
misalignment increases, the spread crosses the equatorial plane of the Bloch sphere. At 50 %, 
the sign of the density shift changes, and therefore the portion of the spread centered around 
this plane averages to zero. The measured density shift is then reduced. 
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Materials and methods 

To first cool and trap the atoms, ^^Sr is loaded into a magneto-optical trap (MOT) from a slowed 
atomic beam. Fig. Sla shows a diagram of relevant energy levels used in cooling, trapping and 
detecting ^^Sr. The MOT is operated on the strong ^Sq— ^Pi transition, and cools the atoms to 
millikelvin temperatures. A second MOT stage operated on the narrow ^Sq— ^Pi transition is 
then used to cool the atoms further to microkelvin temperatures. This second stage MOT uses 
dual-frequency narrow line cooling (i, 2). During the MOT phase the one-dimensional optical 
lattice is superimposed on the atoms in the vertical direction. After the second stage MOT, the 
MOT cooling beams and magnetic gradient are turned off, leaving ~2000 atoms trapped in the 
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optical lattice. After the lasers and magnetic fields are turned off, the temperature of atoms 
remaining in the lattice is ~ 3 fiK. 

To vary the temperature of the atoms, additional cooling or heating is performed on the 
lattice-trapped atoms. The cooling is performed in three dimensions, with Doppler cooling 
applied along the transverse direction of the lattice, and sideband cooling (heating) along the 
longitudinal axis. The laser beams used for cooling (heating) are red (blue) detuned from the 
^5'o(F = 9/2)-^Pi(F = 11/2) transition. F = f+ J is the total angular momentum, with I 
the nuclear spin and J the total electron angular momentum. The cooling (heating) beams are 
appUed first along the transverse direction for 20 ms, followed by 30 ms along the longitudinal 
direction. Simultaneous with the sideband cooling (heating) the atoms are spin-polarized in a 
weak magnetic bias field (~ 30 mG) using an optical pumping beam resonant with the ^So(F 
= 9/2)— ^Pi(F = 9/2) transition. By choosing the correct polarization for the optical pumping 
beam, atoms are polarized into either the m^r = +9/2 (uif = magnetic quantum number) or mp 
= -9/2 ground state. To measure the final temperature of the atoms, two different techniques 
are used. From time-of-flight analysis the temperature of both the longitudinal and transverse 
directions can be directly measured. By examining the spectral components of sideband spec- 
troscopy as shown in Fig. lA, the longitudinal temperature can also be accurately extracted. 
The total atom number is determined by fluorescence measurements. The average density is 
then calculated using the trap volume. The trap volume is determined using our measured trap 
frequencies and temperatures and assumes a Maxwell-Boltzmann distribution. 

Orienting the lattice in the vertical direction breaks the energy degeneracy between lattice 
sites, strongly prohibiting inter-site tunneling (5). This ensures that collisions occur only be- 
tween atoms in the same lattice site. The lattice is aligned with a slight angle with respect to 
gravity, to accommodate the vertical MOT beams. The probe beam used for spectroscopy is 
aligned coUinear with the optical lattice, and the polarization of the beam is also parallel to that 
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of the lattice. To minimize inhomogeneity due to probe beam focusing, the waist of the probe 
beam is ~5 times larger than the waist of the lattice focus. 



Supporting text 

For a given atom in the electronic ground state ^Sq and motional state (rix, Uy, Uz), the excitation 
probability during the spectroscopy pulse is 

where A is the detuning of the clock laser from resonance, t is the length of the spectroscopy 
pulse, and the modified Rabi frequency is given by (4): 

Here VIq is the bare Rabi frequency, Ln^ ^^^ are Laguerre polynomials, and ?7x.j/,2 are the Lamb- 
Dicke parameters for the transverse and longitudinal directions (assuming the projection of the 
misalignment on the transverse plane bisects x and y axes): 



_ sin(A^) / h 



Xl\ 2mvz' ^^^^ 

To calculate the density shift, we return to our two-atom model. As discussed in the main 
text we can explicitly calculate from the anti-symmetrized overlap of the two wavefunc- 
tions. The antisymmetrized wavefunction \^ as) is given by 

\^AS) = ^ (^1^2 - ^2^l) (S4) 

= i= {{a\g) + mMa) + S\e)) - {j\g) + S\e)){a\g) + /3|e))) (S5) 
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and 

G?^{a{t),P{t)n{t),S{t)) = (*A5|*A5) (S6) 

= 1- |a7*+/55*|l (S7) 

Each atom has a slightly different flfi, as given in Eq. S2, and 

. . ^Rlt 

(3{t) = cos ^- + z— sm ^- 

= COS— — sm— — , (S8) 

where Ojiji ^ = \J^i,2^ + A^ is the generalized Rabi frequency. The inhomogeneity is a con- 
sequence of the spread in Rabi frequencies ^nx,ny,n^, due to the excitation of different motional 
states. We can characterize this by defining an average Rabi frequency Q and its RMS spread 
AQ, 

(S9) 

'^x '^y 

AQ'^ = Q^-lf. (SIO) 



Here qx,y,z{'i^x,y,z) is the normalized Boltzmann factor. The inhomogeneity in the Rabi fre- 
quencies (AQ) is affected by T through the Boltzmann factor and by A^ through the probe 
misalignment. To approximate the average density shift for the entire atomic ensemble, we set 
f2i = r2 + and = f2 — in Eq. S8 for our two-atom model. At a time t during the 
spectroscopy pulse, the atoms experience an ensemble-averaged shift: 

Ai/(i) = ?^G'SJ(n+ AQ,n- AQ)(p5-Pe)- (Sll) 



This shift evolves during the spectroscopy pulse, and for the final density shift we time-average 

(2) 

G\2 over the total pulse length t^r as 



1 r 

Al^ave = — / Aiy{t)dt. (S12) 
Jo 
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Figure S 1 : Experimental setup for the Sr optical clock, (a) Relevant energy levels for laser 
cooling and trapping are shown. The first MOT stage is operated on the strong transition at 461 
nm. During this MOT stage atoms that are lost to the and ^Pq states are pumped back to 
the ^So state using two repump beams at 679 nm and 707 nm. The second MOT is operated 
on the narrow 689 nm transition, and the wavelength of the clock laser used for spectroscopy is 
698 nm. After the experiment is performed, atoms in the ^Sq state are detected by measuring 
fluorescence on the strong ^Sq— ^Pi transition. Atoms that have been excited to '^Pq by the 
clock laser pulse are transferred back to the ^Sq using the two repumping beams, and are again 
measured by fluorescence, (b) Two sample spectra. Here, ^^Sr atoms are successively polarized 
into either the mp = 9/2 or -9/2 state. The separation between the two spectroscopy lines is 
~250 Hz. During the actual experiment, we polarize the atoms to each pure state successively, 
and lock near the full-width at half-maximum of the peak (indicated by dots). By moving these 
lock points the desired excitation fraction is selected. 
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